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' Abstract 

A class of dynamical dark energy models, dubbed Spinor Quintom, can be 
constructed by a spinor field ip with a nontraditional potential. We find that, if 
choosing suitable potential, this model is able to allow the equation-of-state to 
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cross the cosmological constant boundary without introducing any ghost fields. 



\& ■ In a further investigation, we show that this model is able to mimic a perfect 

o 

fluid of Chaplygin gas with p = —c/p during the evolution, and also realizes 



the Quintom scenario with its equation-of-state across —1. 



X 

1 Introduction 

The recent data from type la supernovae and cosmic microwave background (CMB) 
radiation and so on [lj [2j [3] have provided strong evidences for a spatially flat and 
accelerated expanding universe at the present time. In the context of Friedmann- 
Robert son- Walker (FRW) cosmology, this acceleration is attributed to the domination 
of a component with negative pressure, called dark energy. So far, The nature of dark 
energy remains a mystery. Theoretically, the simplest candidate for such a component 
is a small positive cosmological constant, but it suffers the difficulties associated with 
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the fine tuning and the coincidence problems. So many physicists are attracted by 
the idea of dynamical dark energy models, such as Quintessence [3], Phantom [5], 
K-essence [6], Quintom [7] and so on (see Refs. [9] for a review). Usually, dark 
energy models are constructed by scalar fields which are able to accommodate a 
rich variety of behaviors phenomenologically. However, there is another possibility 
that the acceleration of the universe is driven by a classical homogeneous spinor 
field. Some earlier studies on applications of spinor fields in cosmology are given in 
Refs. [TOl [TT1 IT2] . In recent years there are many works on studying spinor fields as 
gravitational sources in cosmology, for example: see Refs. p31 [H] for inflation and 
cyclic universe driven by spinor fields; see Refs. [15] for spinor matter in Bianchi 
Type I spacetime; see Refs. [TH [TT] for a dark energy model with spinor matter; and 
so on. 

Although the recent fits to the data in combination of the 3- year WMAP [18] , the 
recently released 182 SNIa Gold sample [19] and also other cosmological observational 
data show remarkably the consistence of the cosmological constant, it is worth noting 
that a class of dynamical models with the equation-of-state (EoS) across —1 Quintom 
is mildly favored [2UJ [2TJ . In the literature there have been a lot of theoretical studies 
of Quintom-like models. For example, motivated from string theory, the authors of 
Ref. [22J realized a Quintom scenario by considering the non-perturbative effects of a 
generalized DBI action. Moreover, a No-Go theorem has been proven to constrain the 
model building of Quintom [23] (see Ref. [23] for some earlier considerations), and 
according to this No-Go theorem there are models which involve higher derivative 
terms for a single scalar field |25j, models with vector field [26], making use of an 
extended theory of gravity [27J, non-local string field theory |28j . and others (see e.g. 
[291 1301 l3T]). The similar work applied in scalar-tensor theory has also been studied 
in Ref. [321 [33]. 

Usually, a Quintom model involves a ghost field with its kinetic term to be negative 
which leads to quantum instability. In this paper we study the dark energy model 
with spinor matter. Interestingly, we find that this type of model can realize the 
Quintom scenario with its EoS across the cosmological constant boundary w = —1 
without introducing a ghost field. Instead, the derivative of its potential with respect 
to the scalar bilinear ipifj, which is defined as the mass term, becomes negative when 
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the spinor field lies in the Phantom-like phase. If this model can realize its EoS across 
— 1 more than one time, the total EoS of the universe can satisfy w > — 1 during the 
whole evolution which is required by the Null Energy Condition (NEC) [34], and we 
expect to treat this process as a phase transition merely existing for a short while. 
Moreover, due to a perfect mathematic property of the spinor field, it is possible to 
combine the Quintom scenario and the picture of Chaplygin gas with the EoS evolving 
from to —1 smoothly in Spinor Quintom. In the literature a dark energy model of 
Chaplygin gas has been proposed to describe a transition from a universe filled with 
dust-like matter to an accelerated expanding stage, and hence it has been argued that 
the coincidence problem of dark energy may be alleviated in this model [35] . 

Interestingly, in our model we are able to evade the drawbacks of considering 
the Phantom field, of which the kinetic energy is negative and so is unstable in 
quantum level. We notice that Ref. [36] has investigated the quantum stability of 
a Phantom phase of cosmic acceleration and shown that a super-acceleration phase 
can be obtained by quantum effects. Our model is different from that one since 
the super-acceleration is realized by the background contribution. However, it is 
interesting that both the two approaches are stable in quantum level since the first 
order of perturbation theory can be defined. 

This paper is organized as follows. In Section II, we simply review the basic algebra 
of a spinor field in FRW universe which is minimally coupled to Einstein's gravity. In 
Section III, we present the condition for the spinor field to realize a Quintom scenario 
and give some detailed examples. In Section IV, we provide a unified model of Spinor 
Quintom and a perfect fluid of Chaplygin gas by taking certain potentials. By solving 
the model numerically, we will study the evolution of its EoS and fraction of energy 
density. Section V is the conclusion and discussions of our paper. 

2 Algebra of A Spinor Field 

To begin with, we simply review the dynamics of a spinor field which is minimally 
coupled to Einstein's gravity (see Refs. [371 ES, [39] for detailed introduction). Fol- 
lowing the general covariance principle, a connection between the metric g^ u and the 
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vierbein is given by 



9liv e a e b ~ Vab 



where denotes the vierbein, g^ v is the space-time metric, and r\ a \> is the Minkowski 
metric with r] a b = diag(l, — 1, — 1, — 1). Note that the Latin indices represents the 
local inertial frame and the Greek indices represents the space-time frame. 
We choose the Dirac-Pauli representation as 

^f 1 ), y=f *'), (2) 

\ o -l y \-<Ti o J 

where <j{ is Pauli matrices. One can see that the 4x47° satisfy the Clifford 
algebra{7 a , r y b } = 2r] ab . The 7° and provide the definition of a new set of Gamma 
matrices 

r" = < 7 Q , (3) 

which satisfy the algebra {T^, T u } = 2g fJ/U . The generators of the Spinor representation 
of the Lorentz group can be written as S afc = ^[7 a ,7 fc ]. So the covariant derivative 
are given by 

= (d^ + n^, (4) 

D^ip = d$ - ipVL^ , (5) 

where the Dirac adjoint ip is defined as The 4x4 matrix = ^u^ a bT l ab is the 

spin connection, where uj^ = e^V are Ricci spin coefficients. 

By the aid of the above algebra we can write down the following Dirac action in 
a curved space-time background [HI [161 STJ] 

= J d*x e [faVD^ - D^Y^) - V] . (6) 

Here, e is the determinant of the vierbein and V stands for the potential of the 
spinor field tp and its adjoint ?/>. Due to the requirement of covariance, the potential V 
only depends on the scalar bilinear ipip and "pseudo-scalar" term -07 s ijj. For simplicity 
we drop the latter term and assume that there is V = V(ipip). 
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Varying the action with respect to the vierbein e£, we obtain the energy-momentum- 
tensor, 

% — — — — 

= -[i/jT u D^ifj + ipT^D u ifj - D^ifjT u ifj - D V ^Y^\ - g^C^ . (7) 

On the other hand, varying the action with respect to the field ■?/>, tp respectively 
yields the equation of motion of the spinor, 

dV - dV 

i^D^--^ = Q, lD ^ + — = V. (8) 

Note that, we use units 8iiG — h — c — 1 and all parameters are normalized by 
M p = l/V8nG in the letter. 



3 A universe driven by Spinor Quintom 

3.1 Dynamics of a spinor field 

In this paper we deal with the homogeneous and isotropic FRW metric, 

ds 2 = dt 2 -a 2 {t)dx 2 , (9) 

where a stands for the scale factor and we choose today's scale factor a — 1. Corre- 
spondingly, the vierbein are given by 

= ef = V (10) 
a 

Assuming the spinor field is space- independent, the equation of motion reads 

3 



ip + ^Hip + i^v'ip = o, (ir 
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i\) + -Hip - i 7 °V r V = , (12) 

where a dot denotes a time derivative 'J^' and a prime denotes a derivative with 
respect to ip-ip, and H is the Hubble parameter. Taking a further derivative, we can 
obtain: 

N 

W = -3 , (13) 
a 



where N is a positive time-independent constant and we define it as today's value of 

From the expression of the energy-momentum tensor in Eq. (jTJ), we get the energy 
density and the pressure of the spinor field: 

= T° = V, (14) 
Pi> = -Tt = V'W>-V, (15) 

where Eqs. ( ITT!) and (|T2|) have been applied. The EoS of the spinor field, defined as 
the ratio of its pressure to energy density, is given by 

w^ = -f = -l + — — . (16) 

Simply taking the potential to be power- law-like V = V (^) a with a as a nonzero 
constant, we obtain a constant EoS: 

w$ — — 1 + a . (17) 

In this case, the spinor matter behaves like a linear-barotropic-like perfect fluid. For 
example: if a = |, we can get p^ ~ a~ 4 and w$ = \, which is the same as radiation; 
if a = 1, then ~ a -3 and = 0, this component behaves like normal matter. 

Furthermore, the spinor matter is able to realize the acceleration of the universe 
if a < |. So it provides us a potential motivation to construct a dynamical dark 
energy model with the spinor matter. Moreover, as introduced at the beginning of 
Section I, there is evidence in the recent observations to mildly support a Quintom 
scenario with the EoS of dark energy across —1. In the following, we emphasize our 
investigation on constructing Quintom dark energy model with the spinor field, which 
is called Spinor Quintom. 



3.2 Evolutions of Spinor Quintom 

To keep the energy density positive, one may see that there is uy, > — 1 when V > 
and Wj/j < — 1 when V < from Eq. (fl6l) . The former corresponds to a Quintessence- 
like phase and the latter stands for a Phantom-like phase. Therefore it requires the 
derivative of the potential V to change its sign if one expects a Quintom picture. In 
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terms of the variations of V, it shows that there exists three categories of evolutions 
in Spinor Quintom: 

(i) , there is 

V > -> V < , 

which gives a Quintom-A scenario by describing the universe evolving from Quintessence- 
like phase with w$ > — 1 to Phantom-like phase with uy, < — 1; 

(ii) , there is 

V < o -> V > o , 

which gives a Quintom-B scenario for which the EoS is arranged to change from below 
— 1 to above — 1; 

(iii) , V changes its sign for more than one time, then one can obtain a new Quintom 
scenario with its EoS crossing —1 many times, dubbed Quintom-C scenario. 

In the following, we will take different potentials of Spinor Quintom to provide 
the three kinds of evolutions mentioned above 0. 

To begin with, we shall investigate Case (i) and provide a Quintom-A model. 
We use the form of potential V = Vo[{jpip — b) 2 + c], where Vq, b, c are undefined 
parameter. Then we get V = 2Vo('0'0 — &) an d the EoS: 

( ^)2 _ b 2 _ c 

* (V>V) 2 - 2fo/>V + b 2 + c' { ' 

According to Eq. ffl3l) . one finds that ifjip is decreasing along with an increasing 

scale factor a during the expansion of the universe. From the formula of V, we 

deduce that at the beginning of the evolution the scale factor a is very small, so ipip 

becomes very large and ensures V > at the beginning. Then tpif) decreases along 

with the expanding of a. At the moment of i/jip = b, one can see that V = which 

results in the EoS = — 1. After that V becomes less than 0, the universe enters 

a Phantom-like phase. Finally the universe approaches a de-Sitter spacetime. This 

-'^Note that we choose the potentials phenomenologically without any constraints from quantum 
field theory or other consensus. From the phenomenological viewpoint, this is okay if we treat the 
background classically while deal with the perturbations in quantum level, just as what is done in 
inflation theory. 
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behavior is also obtained by the numerical calculation and shown in Fig. [TJ From 
this figure, one can read that the EoS w$ starts the evolution from 1, then mildly 
increases to a maximum and then begin to decrease. When ipip = b, it reaches the 
point w^, = — 1 and crosses — 1 from above to below smoothly. After that, the EoS 
sequentially decreases to a minimal value then increases and eventually approaches 
the cosmological constant boundary. 




Figure 1: Plot of the evolution of the EoS in Case (i) as a function of time. In the 
numerical calculation we take the potential of the spinor field as V = Vo[(ipip — b) 2 + c\, 
where we choose Vo = 1.0909 xl0~ 117 , b = 0.05 and c = 10~ 3 for the model parameters. 
For the initial condition we take N = 0.051. 



In Case (ii) we choose the potential as V = Vo[— (ipip — b)ipip + c]. Then one can 
obtain V = V (-2ipip + b) and the EoS 

= -WO 2 - c 
Wf -(ipip) 2 + bipip + c' 1 ' 

Initially V is negative due to the large values of ipip. Then it increases to when 
ipip = | whereafter changes its sign and becomes larger than 0, in correspondence 
with the Case (ii), dubbed Quintom-B model. From Fig. [2j we can see that the EoS 
evolves from below —1 to above —1 then finally approaches to —1. 

In the third case we explore a Quintom scenario which gives the EoS across — 1 
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Figure 2: Plot of the evolution of the EoS in Case (ii) as a function of time. In the 
numerical calculation we take the potential of the spinor field as V = Vo[—(ipip — 
btyip + c], where we choose Vq = 1.0909 x 10~ 117 , b = 0.05 and c = 10~ 3 for the model 
parameters. For the initial condition we take N = 0.051. 



for two times in virtue of the potential 

V = V [$ip-b) 2 <ipip + c] . (20) 
From the expression of the potential, one have 

V = V (^-b)(3^-b) , (21) 

and the EoS 

= 2(^) 3 -26(#) 2 -c 
W ^ (^) 3 -26(^) 2 + & 2 (^>) +c ' 
Evidently the equation V — has two solutions which are ipip = b and iptp = |, thus 
V changes its sign two times. From the expression of the EoS, we find that > — 1 
in the beginning. When the value of ipip equals to b, it crosses —1 for the first time. 
After the first crossing, it enters the phantom-like state and continuously descends 
until passes through its minimum, then ascends to ipip = | and then experiences the 
second crossing, and eventually moving up to the Quintenssence-like phase. This is 
shown in Fig. [31 One can see that the big rip can be avoided in this case. 
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Moreover, taking the component of normal matter with the energy density p m oc 
1/a 3 into consideration and the EoS w m = 0, we can see that the EoS of the universe 
w u = ^ p satisfies the relation w u > — 1 in this case. As is argued in Ref. [31], 
NEC might be satisfied in the models though < — 1 only stays for a short period 
during the evolution of the universe. 




Figure 3: Plot of the evolution of the EoS in Case (iii) as a function of time. The 
magenta solid line stands for the EoS of Spinor Quintom and the green dot line stands 
for that of the whole universe. In the numerical calculation we take the potential of 
the spinor field as V = Vo[(^V ; — b) 2 %pip + c], where we choose Vo = 1.0909 x 10~ 117 , 
b = 0.05 and c = 10~ 3 for the model parameters. For the initial condition we take 
N = 0.051. 



One important issue of a dark energy model is the analysis of its perturbations. 
To study this issue, we might be able to learn to what degree the system is stable 
both in quantum and classical level. Usually systems with w < — 1 show some nasty 
instabilities, for example, a Phantom universe suffers a Big rip singularity. Although 
this singularity can be avoided in most Quintom models which usually enter a de- 
Sitter expansion in the final epoch, all the scalar models of Quintom by now suffer 
a quantum instability since there are negative kinetic modes from ghost fields. Here 
we would like to show the perturbation theory of Spinor Quintom crudely. Since we 
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do not introduce any ghost fields in our model, the crossing of phantom divide is 
achieved by the spinor field itself and does not perform any particular instabilities. 
In order to simplify the derivative, we would like to redefine the spinor as ipN = 

3 

a 2-0- Then perturbing the spinor field, one gives the perturbation equation as follows, 

d 2 

-T^SipN - V 2 6ip N + a 2 [V 2 + ij°{HV - 3HV"W)] 5?p N 

(Jj I 

= -2a 2 V'V"5{i>^ N - Wdr[aV"8(W>)\*l>N , (23) 

where r is the conformal time defined by dr = dt/a. Since the right hand side of the 
equation decays proportional to a~ 3 or even faster, we can neglect those terms during 
the late time evolution of the universe for simplicity. 

From the perturbation equation above, we can read that the sound speed is equal 
to 1 which eliminates the instability of the system in short wavelength. Moreover, 
when the EoS w crosses —1, we have V = at that moment and the eigen function 
of the solution to Eq. (1231) in momentum space is a Hankel function with an index 
|. Therefore, the perturbations of the spinor field oscillate inside the hubble radius. 
This is an interesting result, because in this way we might be able to establish the 
quantum theory of the spinor perturbations, just as what is done in inflation theory. 

Note that, the above derivative does not mean that Spinor Quintom is able to 
avoid any instabilities. We still do not study the effects of the right hand side of 
Eq. ( |23i) which may destroy the system under some certain occasions. Another 
possible instability may be from the quantum effect that our model is unable to be 
renormalized. The more detailed calculation will be investigated in the future works. 

4 A unified model of Quintom and Chaplygin gas 

In the above analysis, we have learned that a spinor field with a power-law-like po- 
tential behaves like a perfect fluid with a constant EoS. However, it is still obscure to 
establish a concrete model to explain how a universe dominated by matter evolves to 
the current stage that is dominated by dark energy. In recent years, another interest- 
ing perfect fluid with an exotic EoS p — —c/p has been applied into cosmology [35] in 
the aim of unifying a matter dominated phase where p oc 1/a 3 and a de-Sitter phase 
where p = —p which describes the transition from a universe filled with dust-like mat- 
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ter to an exponentially expanding universe. This so-called Chaplygin gas [35] and its 
generalization [31] has been intensively studied in a literature. Some possibilities for 
this model motivated by field theory are investigated in [42J. A model of Chaplygin 
gas can be viewed as an effective fluid associated with D-branes [Ml EH], an d also can 
be obtained from the Born-Infeld action [321 HE] . The combination of Quintom and 
Chaplygin gas has been realized by the interacting Chaplygin gas model [37J as well 
as in the framework of Randall-Sundrum braneworld [48J . 

The Chaplygin gas model has been thoroughly investigated for its impact on the 
cosmic expansion history. A considerable range of models was found to be consistent 
with SN la data [31], the CMBR [50], the gamma-ray bursts [51], the X-ray gas mass 
fraction of clusters [52J, the large scale structure [53], and so on. 

Here, we propose a new model constructed by Spinor Quintom which combines 
the property of a Chaplygin gas. The generic expression of the potential is given by 



V 



+ c , (24) 



where f{ipip) is an arbitrary function of Altering the form of f(ipip), one can 
realize both the Chaplygin gas and Quintom scenario in a spinor field. 

Firstly, let us see how this model recovers a picture of generalized Chaplygin gas. 
We take /(V>V) = ^(ipip)^ 13 , an d then the potential is given by 



V= ^V (^y+P + c. (25) 
Due to this, we obtain its energy density and pressure 



= 1+ Vv (wy+e + c, (26) 



pip = -c\y (w>y+ f, + c]-"z . (27) 

Now it behaves like a generalized Chaplygin fluid which satisfies the exotic relation 
Pip — — -T- To be more explicitly, we take (3=1, then get the expressions of energy 
density and pressure to be 



N 2 c 
P, = ^ + c, * = (28) 

In this case a perfect fluid of Chaplygin gas is given by a spinor field. Based on 
the above analysis, we may conclude that this simple and elegant model is able to 
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mimic different behaviors of a perfect fluid and so that accommodate a large variety 
of evolutions phenomenologically. 

In succession, we will use this model to realize a combination of a Chaplygin gas 
and a Quintom-A model which is mildly favored by observations. Choosing f(ipip) to 
be f(ipip) = Voi^pip ~ b) 2 , we get the potential 

V = yV (^-&) 2 + c , (29) 

where Vo, b, c are undetermined parameters. So we obtain the derivative of the 
potential 

- (30) 
s/V (W ~ W 2 + c 

and the EoS 

m * = " 1 + y.W-t)' + c ' (31) 

respectively, and the crossing over —1 takes place when = 6. 

During the expansion is decreasing following Eq. ({TBI) . From the formula 
of the EoS 0311) . we deduce that at the beginning of the evolution the scale factor 
a — > so ipip — > oo. To neglect the terms of lower order of and the constants, 
the EoS at early times evolves from which describes the era of matter dominated. 
Along with the evolution, uy, increases from to the maximum and then starts to 
decrease. When ipif) = b the EoS arrives at the cosmological constant boundary 

= — 1 and then crosses it. Due to the existence of c term, the EoS approaches the 
cosmological constant boundary eventually. In this case the universe finally becomes 
a de-Sitter space-time. Considering a universe filled with normal matter and such a 
Spinor Quintom matter, we take the numerical calculation and show the evolution of 
the EoS in Fig. HI Moreover, we display the evolvement of fraction densities of the 
normal matter Q m = p m /(Pm + Pip) and Spinor Quintom = p^/(p m + p^) in Fig. [5j 
It is evident that there is an exact ratio of these two components from the beginning 
of evolution, in relief of fine-tuning problems and accounting for the coincidence 
problem. Then they evolve to be equal to each other in late times. Along with the 
expansion of a, the dark energy density overtakes the matter energy density driving 
the universe into an accelerating expansion at present and eventually dominates the 
universe completely, which describes an asymptotic de-Sitter spacetime. 
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Figure 4: Plot of the EoS of the unified model in Eq. (1291) as a function of time. In 
the numerical calculation we take V = 3.0909 x 10 -239 , b = 0.05 and c = 9 x 10~ 241 . 
For the initial conditions we take N = 0.051. 



Q 




Figure 5: Plot of the fraction of energy density of dark energy (the violet solid line) 
and normal matter (the red dash line) as a function of time. 
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5 Conclusion and Discussions 



The current cosmological observations indicate the possibility that the acceleration 
of the universe is driven by dark energy with EoS across —1, which will challenge 
the theoretical model building of the dark energy if confirmed further in the future. 
In this paper we have studied various Quintom scenarios in virtue of a spinor field 
and proposed a unified model of Spinor Quintom and a generalized Chaplygin gas. 
As shown in the present work, this model can give rise to the EoS crossing the 
cosmological constant boundary during the evolution by varying the sign of the term 
V. Compared with other models with w across —1 in the literature, so far the present 
one is also economical in the sense that it merely involves a single spinor field. 
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